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Abstract — In the recent years , the contribution of the 
wind power to energy supply has increased considerably; 
hence, the wind farms have to be able to participate to the 
grid power stability. In this paper, an optimization 
algorithm allows obtaining the reactive and active power 
dispatch in a wind power plant is presented. The aim of 
the proposed algorithm is to minimize the power losses 
and the difference between the reactive power obtained 
and required by the transmission system operator at the 
point of common coupling. The simulation results show 
the validity and the performance of the proposed 
algorithm. 
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I. INTRODUCTION 

In the last years, the integration of wind energy into 
power grids has grown significantly. In order to ensure 
the power quality of the grid, transmission system 
operators (TSO) in different countries require from wind 
farms (WFs) to be able to contribute to ancillary services, 
especially in reactive power control. 

Many works have been done in order to get an optimal 
reactive power dispatch using several methods such as 
particle swarm optimization (PSO) technique. In [1] a 
PSO technique for a reactive power wind farm function is 
presented. This technique allows obtaining the reactive 
power reference for each wind turbine with the aim to 
optimize the reactive power dispatch at a wind farm and 
minimize its power losses. 

This paper is organized as follows: First, the doubly fed 
induction generator (DFIG) is presented. Second, the 
reactive power limitations of a wind turbine are 
determined using the method proposed in [2]. Then, a 
multi objective function that allows getting under several 
constraints both the reactive and active power set point 
for each wind turbine (WT) is presented. Finally, 
simulation results are reported. 


II. DOUBLY FED INDUCTION GENERATOR 

(DFIG) 

The doubly fed induction generator is coupled to the grid 
with the stator windings, while the rotor windings are 
connected to the grid via a back-to-back converter [3], as 
shown in figure 1. The DFIG exchanges the power with 
the grid through the stator windings as well as the rotor 
windings. The main part of the power passes from the 
generator through the stator into the grid, whereas only a 
fraction of the power is passed from the rotor windings 
through the power converter [2] . 

The DFIG technology is the most widely used generator 
in the wind farms for several reasons. First, it has the 
ability to control electrical torque (hence active power) 
and reactive power exchange with the grid. Besides, the 
DFIG is the cheapest solution to realize variable speed 
operation because the converters are sized only for 20% - 
35% of the stator power (not total turbine power) 
depending on the slip (or operating speed) range and 

reactive power requirements [4]. 

• Pr rotor power • • ‘ *l«p 

• Pi : stator power • Pmcch mechanic*) power 



Fig. 1: The basic layout of a DFIG wind turbine [5] 


III. DFIG CAPABILITY LIMITS CURVE 

We use the method proposed in [2] in order to get the 
reactive power capability of a 2MW DFIG based wind 
turbine. In this method, we consider that the reactive 
power capability is limited by three parameters: stator 
current (Is), rotor current (Ir), and rotor voltage (Vr). In 
order to obtain the PQ diagram, the stator voltage is 
considered to be equal to 1 p.u, and the steady state T- 
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IV. PARTICLE SWARM OPTIMIZATION 
METHOD 


equivalent circuit, as shown in figure 2, is used to derive 
the complex powers from the stator and the rotor 
windings. 


Rs X s X r 

s 



Fig. 2: Steady state T-equivalent circuit for the DFIG [2] 

The reactive power capability of a DFIG is obtained by 
the most restrictive of the three limitations. Figure 3 
shows the PQ curve of 2MW DFIG based wind turbine 
using the parameters illustrated in table I. 



11 12 13 M IS 


Particle swarm optimization is a heuristic global 
optimization method put forward originally by Doctor 
Kennedy and E berhart in 1995 (Kennedy J,Eberhart 
R,1 995 ;EberhartR, Kennedy J,1995). It is developed from 
swarm intelligence and is based on the research of bird 
and fish flock movement behavior [7]. 

The particle swarm model consists of a group of particles 
that are randomly initialized in the d-dimensional search 
space. During an iterative process, particles explore this 
space effectively by exchanging information to find the 
optimal solution. Each i-th particle is described by its 
position xi, velocity vi, and best position pbesti. 
Moreover, the particles have access to the best global 
position gbest that has been found by any particle in the 
swarm [8]. 

Then, each particle updates its coordinates based on its 
own best search experience pbesti and gbest according to 
the following velocity and position update equations: [9] 

v k+1 = w vj 5 + cl rl(pbest k - x} 5 ) 

+ c2 r2(gbest k — x k ) 

x k+1 = x k + v k+1 (2) 

Where: 

w: Inertia weight 

cl, c2: Acceleration coefficients 

rl,r2: Two separately generated uniformly distributed 
random numbers in the range [0,1] added in the model to 
introduce stochastic nature. 

The inertia weighting factor for the velocity of particle is 
defined by the inertial weight approach 




Fig. 3: Results reactive power capability used in this 
paper 


Table I: REpower MM82 parameters [6] 


Machine parameter 

Value 

Nominal active power 

2MW 

Nominal stator voltage 

690 V 

Stator resistance 

0.00206 Q 

Stator inductance 

0.032 Q 

Rotor resistance 

0.0028 Q 

Rotor inductance 

0.021 D 

Magnetizing resistance 

36.4 G 

Magnetizing inductance 

0.83 Q 

Turn ratio 

2.43 

slip 

-0.2 


W k = W max - ■ 


x k 


(3) 


Where: 

k max is the maximum number of iterations, and k is the 
current number of iterations. 

w max 2in ^ w min are the upper and lower limits of the 
inertia weighting factor, respectively. 

V. PROBLEM FORMULATION 

We propose to use a multiobjective optimization function 
to calculate the reactive and active power set point for 
each wind turbine within the wind power plant. The 
multiobjective function is expressed as follows: 

min F(X) = | Qp^ - Q pcc \ + X P loss (4) 

Where: 


• X=( Q1,Q2,Q3,Q4,Q5,P1,P2,P3,P4,P5) 


www.ijaers.com 


Page | 16 




International Journal of Advanced Engineering Research and Science (IJAERS) [Vo 1-4, Issue-3, Mar- 2017] 

https://dx.doi.Org/10.22161/iiaers.4.3.2 ISSN: 2349-6495(P) / 2456-1908(0) 


Pi is the WT active power generation. 

Qi is the WT reactive power consumption/generation. 

• Qpcc is the reactive power required at the PCC by the 
TSO 

• Qpcc is the reactive power generated by the wind 
farm at the PCC ,and it is obtained as follows: 

Qpcc — 7 Qgt Qlosscs (^) 


Where Q g . is the generated or absorbed reactive power of 
each iDFIG, and Q LoS ses is the reactive power losses 
within the wind power plant. 

• Pioss is the real power losses within the wind power 
plant. 

• A, is the weight coefficient. 

The minimization of the multiobjective function is subject 
to the following constraints: 

1). The node power equation 


Nb 


p i= Ui ^ Uj ( Gij cos d t j + Bij sin d i} ) 

(6) 

i= i 


Nb 


Qi= Ui^ Uj ( Gij sin Qjj - B i} cos 6 tj ) 

(7) 

J = i 


Where: 



Ui and Uj are the voltage amplitude of nodes i and j 
respectively; 

0ij = 0i - 0j is voltage phase angle difference of node i and 

j; 

Gij and Bij are mutual conductance and susceptance of 
admittance matrix respectively; 

Pi and Qi are injected active and reactive power of node i. 

2) . DFIG reactive capability limits 

Qmin < Q. < Qmax ( 8 ) 

3) . DFIG active power limits 

pmin < p < pmax (9) 

The PSO algorithm is initialized with the population of 
individuals being randomly placed between the space of 
possible values [Q^,Q^ X ] and [P^,P ^ X ], and it 
looks for the optimal solution by updating individual 
generations. The velocity and the position of each particle 
are updating, at each iteration, according to its previous 
best position (Pbest) and the best of the group gbest, as 
illustrated in figure 4. 



Fig. 4 : Flow chart of the proposed PSO algorithm 


VI. SIMULATION RESULTS 

An example system is used as shown in figure5 in order 
to verify the validity and performance of the proposed 
optimization algorithm. The wind farm consists of five 2 
MW wind turbines, and it is connected to a 20 kV 
distribution system which exports power to 63 kV grid 
through a 10 km 20 kV feeder. 



Fig. 5: An example system for simulation 
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Three different cases to do the reactive power dispatch 
are tested and compared. The WTs within the wind power 
plant are operating at full active power in both case 1 and 
2; also, the reactive power Q required at the PCC by the 
TSO is equal to 0,7 Mvar in the all cases. 

• Casel: the reactive power required at the PCC is 
proportionally distributed between the WTs within 
the wind power plant. 

• Case2: The reactive power reference for each WT is 
obtained using the PSO optimization technique 
proposed in [1], where the OF is to minimize both the 
deviation between the reactive power Q^ e f c required 
atthe PCC and the reactive power generated by the 
WF and the losses along the branches of the WF, as 
shown in equation (4). In this case X= (Ql, Q2, Q3, 
Q4, Q5). 

Three values of the weight coefficient A are tested: 

X = 0 ;X = 0,5 ; X = 0,9 

• Case3: the reactive and active power set point for 

each wind turbine within the wind power plant are 
obtained using the objective function as shown in 
equation (4). In this case 

X=(Q1,Q2,Q3,Q4,Q5,P1,P2,P3,P4,P5). 

The same values of A,used in the case 2 are tasted in 
this case. 

The results of the reactive and active power dispatching 
in the three cases are presented in table II: 


Table II: Simulation results 



Cas 


Case2 



Case3 



el 

x=o 

X=0, 

5 

x=o, 

9 

A,=0 

X=0, 

5 

x=o, 

9 

Ql 

0,1 

0,46 

0,46 

0,46 

0,13 

0,38 

0,17 

Q2 

4 

0,46 

0,46 

0,46 

- 

0,04 

0,33 

Q3 

0,1 

0,46 

0,46 

0,46 

0,01 

0,07 

0,21 

Q4 

4 

0,46 

0,46 

0,46 

0,29 

0,28 

0,23 

Q5 

0,1 

0,46 

0,46 

0,46 

0,20 

0,12 

- 

PI 

4 

2,00 

2,00 

2,00 

0,28 

1,00 

0,04 

P2 

0,1 

2,00 

2,00 

2,00 

1,00 

1,00 

1,00 

P3 

4 

2,00 

2,00 

2,00 

1,00 

1,00 

1,00 

P4 

0,1 

2,00 

2,00 

2,00 

1,00 

1,00 

1,00 

P5 

4 

2,0 

0 

2,0 

0 

2,0 

0 

2,0 

0 

2,0 

0 

2,00 

2,00 

2,00 

1,00 

1,00 

1,00 

1,00 

1,00 

Object 

- 

0,01 

0,09 

0,16 

0,00 

0,00 

0,01 
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ive 

functi 

on 


82 

82 

23 

44 

98 

77 

| o set 

I Ypcc 

— Qpcc 1 

1,5 

9 

0,01 

82 

0,01 

82 

0,01 

82 

0,00 

44 

3,41 

e-7 

3,12 

e-7 

Line 

0,1 

0,16 

0,16 

0,16 

0,01 

0,01 

0,01 

losses 

(MW) 

61 

01 

01 

01 

96 

96 

96 


The case 1 has the most important value of line losses, but 
it decreases significantly in the case 3 which represents a 
reduction of 87,83% in comparison with the case 1. 

We obtain the important value of the error in reactive 
power at the PCC in the case 1. This error decreases in 
the case 2 which represents a reduction of 98,85% in 
comparison with the case 1 . However, the error is almost 
equal to 0 in the case 3 for all the values of A,. 


0,2 . 0.161 0,1601 



Cl C2 C3 
X X 

Fig. 6: Power losses in the WF in each case for 1=0 
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Fig. 7: Power losses in the WF in each case for 7=0,5 
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Fig. 8: Power losses in the WF in each case for 1=0,9 



Cl Cl Ci 


Fig. 9: Error in reactive power at the PCC in each case 
for X=0 
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Fig. 10: Error in reactive power at the PCC in each case 
for X=0, 5 
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Fig. 11: Error in reactive power at the PCC in each case 
for X=0, 9 

Comparing the results of the three cases, we can conclude 
that the case 3 is the most advantageous because it allows 


reducing significantly the line losses in the WF, and 
obtaining an accurate value of the reactive power 
generated at the PCC that corresponds to required by 
the TSO, in comparison with all the other cases. 

VII. CONCLUSION 

In this work, we present an optimization algorithm of 
reactive and active power dispatch in a wind farm. The 
proposed algorithm allows minimizing both the power 
losses in the WF and the deviation between the reactive 
power reference the reactive power generated by 

the wind power plantat the point of common 
coupling Q pcc . By analyzing the obtained results, we can 
conclude that the proposed method is the most suitable to 
use in order to get an accurate value of the reactive power 
at the point of common coupling and minimize the power 
losses in the wind farm. 
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